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SUMMARY 


The  electric  conductivity  tensor  of  a beam-plasma  interacting 
system  is  an  extremely  significant  operator.  Within  the  model  we 
have  chosen,  one  calculates  the  effects  of  the  beam  on  the  plasma 
as  corrections  to  the  classical  results  found  from  Drude  theory. 
This  work  has  been  supported  by  the  Independent  Research  Program 
of  the  Center. 
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I.  INTRODUCTION 

Many  approaches  to  the  understanding  of  the  frequency- 
dependent,  electrical  conductivity  tensor  of  a plasma  have  been 
made  over  the  years.  Depending  on  the  plasma  model  one  was  con- 
sidering, investigators  have  chosen  the  relevant  transport  equa- 
tion that  best  suited  their  needs.  One  of  the  earliest  treat- 
ments of  the  conductivity  was  that  done  by  Drude.^  He  considered 
a gas  of  electrons  that  were  free  to  move  and  collide  with  a 
collection  of  fixed  ions  at  rest.  When  an  external  electric 
field  is  applied  the  motion  of  individual  electrons  is  governed  by 
Newton's  law  which  contains  a field-dependent  term  and  a collision 
frequency  term.  The  resulting  conductivity  depends  inversely  on 
the  collision  frequency  in  the  low  frequency  limit  of  the  applied 
electric  field.  This  is  the  so-called  DC  conductivity.  Since 
the  average  collision  frequency  is  proportional  to  the  mean  velo- 
city of  the  electron,  the  conductivity  goes  inversely  as  the 

square  root  of  the  absolute  temperature  of  the  electrons.  Cohen, 

2 

Spitzer  and  Routley  used  the  Boltzmann  equation  to  obtain  the 
velocity  distribution  function  for  an  electron  gas  mixed  with 
singly-ionized  atoms.  They  treated  the  collision  term  (i.e.  the 
time  rate  of  change  in  the  distribution  function  due  to  colli- 
sional  processes)  as  composed  of  two  parts.  The  first  part  is 

3 

the  traditional  representation  for  close  encounters  which  pro- 
duce large  deflections  in  particle  trajectories.  The  second 
represents  a diffusion  process,  such  as  one  finds  in  Brownian 
motion.  Thus  distant  encounters  produce  only  minor  deflections 
in  the  motion  of  particles.  The  DC  conductivity  obtained  from 
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this  theory  exhibits  a temperature  behavior  that  is  proportional 

4 

to  the  3/2  power.  In  a later  paper,  Spitzer  and  Harm  included 

the  effects  of  electron-electron  interactions  on  the  collision 

term  represented  by  the  Fokker-Planck  equation,  for  a fully 

ionized  gas.  The  results  were  quantitatively  comparable  to  those 

3 5 6 

of  Chapman  and  Cowling  and  Cowling.  Jackson  was  able  to  com- 
bine the  Drude  theory  with  an  expression  for  the  mean  ion-electron 
collision  frequency.  This  expression  results  from  a careful 
treatment  of  the  multiple  scattering  of  electrons  as  they  pass 
through  a finite  region  of  matter.  His  results  compare  with 
those  of  Spitzer  to  within  a factor  of  2.  Paramount  to  this  re- 
sult was  the  assumption  of  small-angle  Coulomb  scattering  as  the 
dominant  process.  These  results  are,  of  course,  for  the  zero 
frequency  limit  (or  DC)  conductivity.  Rosenbluth,  MacDonald  and 

n 

Judd  found  a series  representation  for  the  distribution  function 
satisfying  the  Fokker-Planck  equation.  The  series  is  given  in 
terms  of  Legendre  functions  relative  to  some  axis  of  symmetry 
that  is  specified  by  an  external  field.  The  first  term  in  the 

O 

series  corresponds  to  Chandrashekhar's°  calculation,  the  second, 

2 

to  that  of  Cohen,  Spitzer  and  Routly.  Retaining  higher  order 
terms  in  the  distribution  function  expansion  will  lead  to  a 
temperature  dependence  in  the  DC  electric  conductivity  that  de- 
parts from  the  3/2  power  behavior.  These  approaches  to  the 
electrical  conductivity  do  not  retain  the  dynamics  of  the  trans- 
port process  itself  except  for  the  very  simplistic  Drude  theory. 
The  response  to  some  time-dependent,  external  electric  field,  in 
general,  will  depend  on  the  history  of  the  disturbance  as  well  as 


4 


r 


NSWC/WOL  TR  77-146 

its  present  value.  For  a plasma  medium  that  is  subjected  to  a 

time-varying  external  electric  field,  the  medium  takes  a finite 

amount  of  time  to  respond  to  the  field  and  cannot  keep  up  with 

its  changes.  Thus  the  electric  polarization  of  the  medium  at  any 

given  time  will  be  responding  to  excitations  experienced  at  some 

time  in  the  past.  More  concisely,  the  electric  displacement  and 

the  electric  field  will  be  functionally  related  through  a con- 

9 10 

volution  in  time.  ' Thus  the  electric  conductivity  and  the 
electric  permeability  of  the  medium  are  frequency-dependent 
quantities . 

An  alternative  approach  to  the  study  of  such  time-dependent 
phenomena  is  a theory  of  irreversible  processes  for  a system 
that  is  not  very  far  from  thermodynamic  equilibrium.  This 
theory  relates  the  thermodynamic  average  of  the  time  fluctuations 
of  generalized  forces  which  act  on  a system  to  some  dissipation 
' parameter  or  kinetic  coefficient.  For  the  situation  at  hand,  the 

■ generalized  force  would  be  the  applied,  time-dependent  electric 

field,  the  kinetic  coefficient  would  be  given  by  the  electric 
conductivity.  Nyquist^^  introduced  the  notion  of  coupling  be- 
tween the  voltage  fluctuations  in  a linear  electrical  network  to  the 
electrical  resistance  of  the  system.  The  dissipation  process  is 
the  result  of  the  system  attempting  to  come  into  thermal  equili- 
brium with  the  external  source.  The  fluctuations  in  the  external 
source  serve  as  the  mechanism  for  the  transfer  of  energy  to  the 
dissipative,  or  electrical  system,  and  dissipation  results  from 
the  randomizing  effect  of  the  fluctuations  on  the  external  source. 
Kirkwood  found  an  expression  for  the  friction  constant  that 
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plays  an  important  role  in  the  theory  of  Brownian  motion  through 
the  Langevin  equation.  The  correlation  of  intermolecular  forces 
acting  on  a given  particle  in  some  system  of  particles  is  propor- 
tional to  the  friction  constant.  Callen  and  Welton^^  generalized 
Nyquist's  results  and  obtained  a relationship  between  the  general 
impedance  of  a linear  system  and  fluctuations  of  generalized 
forces.  The  time  correlations  between  a generalized  force  and 

any  function  of  momenta  and  coordinates  for  a statistical  system 

14 

was  considered  by  Takahasi.  His  results  served  to  unify  the 
static  statistical  theory  of  Gibbs^^  to  the  time  dependent 
phenomena  represented  by  thermal  noise  in  electrical  circuits, as 

, discussed  by  Nyquist.^^ 

ii  < 

[ The  first  attempt  to  apply  the  general  concept  of  the  flue- 

I 

tuation  dissipation  theorem  to  the  problem  of  electric  conduction 

in  metals  was  made  by  Kubo  and  co-workers . They  obtained 

^ expressions  for  the  generalized  susceptibilities  for  systems  that 

i are  strongly  coupled  to  external  generalized  forces,  and  related 

these  susceptibilities  to  fluctuations  in  observable  macroscopic 

quantities  such  as  particle  coordinates  or  momenta.  In  the  case 

of  the  conduction  problem,  the  applied  electric  field  plays  the 

role  of  the  generalized  force.  H.  Mori^^  successfully  unified 

12 

Kirkwood's  theory  of  statistical  mechanics  with  the  stochastic 

O 

theory  of  Brownian  motion  embodied  by  the  Langevin  equation. 

The  resulting  generalized  Langevin  equation  successfully  treats 
the  behavior  of  transport  processes  through  the  time  evolution 
of  microscopic  internal  processes  such  as  the  temperature  gradient 
and  stress  tensor  in  fluids.  He  also  distinguishes  between  two 
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types  of  relaxation  processes  in  fluids,  a macroscopic  process 

which  is  represented  by  hydrodynamic  equations  and  a microscopic 

process  which  allows  for  local  thermodynamic  equilibrum.  The 

1 6 

electric  conduction  problem  studied  by  Kubo  is  analogous  to  the 

1 7 

macroscopic  process  in  fluids  studied  by  Mori.  A study  of 

plasma  transport  coefficients  using  the  correlation  function 

1 8 

approach  of  Kubo  and  Mori  was  made  by  Bartis  and  Oppenheim. 

In  particular,  the  self-diffusion  coefficient  was  studied  in  the 
absence  of  external  electric  or  magnetic  fields.  Short  range 
forces,  as  represented  by  a binary  collision  operation,  and  long 
range  Coulombic  forces  are  used  to  represent  the  forces  acting  on 
the  particles  in  the  system.  Diagrammatic  techniques  are  em- 

I 

ployed  to  conveniently  sum  the  non-zero  contributions  to  the 

interactions  and  to  arrive  at  convergent  values. 

Questions  have  been  raised  concerning  the  failure  of  Kubo's 

^ theory  to  distinguish  between  the  observed,  macroscopic  field, 

19 

i and  the  applied  field.  Izuyama  included  the  effects  of  Coulomb 

interaction  between  the  electrons  in  the  metal  as  well  as  the 
interaction  of  the  electrons  with  the  external  field.  His 
spatially-dependent , and  frequency-dependent  conductivity  has  as 
its  leading  term,  the  Kubo  result. Higher  order  terms  corre- 
spond to  the  Coulomb  interaction.  He  obtains  an  expression  for 
the  conductivity  that  represents  the  ratio  of  the  current  res- 
ponding to  the  effective  field  in  the  metal  to  the  macroscopic, 
effective  field.  This  effective  field  is  nothing  more  than  the 
applied  field  plus  the  induced  polarization  field. 

We  shall  use  the  correlation  function  method  developed  by 
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Kubo^^  to  Study  the  electric  conductivity  of  a plasma  in  an 

external  field.  The  origin  of  the  field  will  be  a line  source  of 

charge  and  current,  which  serves  to  model  an  electron 

beam.  The  time  evolution  operator  for  the  current  flowing  in 

the  plasma  shall  be  represented  by  the  usual  one  established  by 
12 

Kirkwood,  the  harmonic, exponential,  Liouville  operator.  A mean 
field  theory  approximation  will  be  introduced  to  simplify  the 
expressions  for  the  interactions  between  the  various  charged 
particles  in  the  plasma.  The  time  dependent  response  function 
is  expressed  in  terms  of  the  radial  and  axial  current  components 
within  the  plasma.  The  Courier  transform  of  this  quantity  directly 
yields  the  tensor  components  of  the  frequency  dependent  electric 
conductivity.  Ignoring  spatial  changes  in  the  Hamiltonian  yields 
the  usual  results  of  Drude  theory,  retention  of  the  changes  in- 
troduces some  additional  complexity  in  the  frequency  dependence 
of  the  conductivity. 

We  shall  postpone  for  a later  paper  the  effects  of  binary 
collisions  on  the  conductivity  tensor,  where  we  shall  study  the 
"nonlocal"  characteristics  of  the  conductivity  as  well. 
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1 

II.  THE  ELECTRIC  CONDUCTIVITY,  THE  ADMITTANCE  FUNCTION 
AND  THE  RESPONSE  FUNCTION 

The  introduction  of  a relativistic  electron  beam  into  a 
neutral  gas  disturbs  the  initial  state  of  the  system  quite  abrupt- 
ly and  dramatically.  If,  however,  instead  of  a sharp  "switching 
on"  of  the  beam,  we  adiabatically  insert  the  beam  into  the  gas 
before  we  define  the  initial  state  of  the  system  and  the  beam 
particle  density  is  much  less  than  the  gas  density,  then  the 
beam-gas  interaction  may  be  represented  by  some  perturbation 
energy.  A simplified  model  of  this  system  shall  be  a plasma  (the 
result  of  ionization  by  the  beam  front)  in  the  presence  of  an  ex- 
ternal current,  charge  source  of  infinite  length  which  is  the 
initial  state  of  the  system  at  t=  - “.  The  details  of  the 
ionization  process  and  other  attendant  phenomena  such  as  recombi- 
nation, excitation,  and  radiative  processes  shall  be  totally  ig- 
nored. What  shall  be  addressed  in  this  paper  is  the  response  of 
the  plasma  to  the  external  field.  A brief  discussion  relating 
the  generated  conductivity  to  the  general  admittance  and  response 
function  will  help  to  establish  the  approach  to  charge  conduction. 

If  we  represent  the  induced  time-dependent  current  as  a 
change  in  some  initial  value  we  can  write 

AJ(t)  = (2tt)  =*  / (2.1) 

» 00 

We  have  not  allowed  for  spatial  variation  in  the  current.  This 
assumption  is  reasonable  if  we  assert  that  the  external  electro- 
magnetic field  varies  quite  slowly  through  space,  over  distances 
large  compared  to  the  Debye  length  of  the  plasma.  Furthermore, 
for  simplicity  we  restrict  ourselves  to  one  dimension  so  that  we 
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may  most  easily  establish  relations  that  can  be  directly  extended 
to  the  three  dimensional  case.  For  convenience,  we  also  assume 
there  is  one  Fourier  component  that  dominates  heavily  over  all 
others,  then  we  may  write  approximately, 

E(t)  ~ (2.2) 

where  E(t)  is  the  external,  perturbing  electric  field. 

We  use  the  generalization  of  Ohm's  law  which  is  valid  for 
rapidly  time  varying  fields  whose  principle  frequencies  are  not 
small  in  comparison  to  those  frequencies  that  characterize  the 
electric  properties  of  the  plasma,  such  as  the  plasma  frequency. 
The  generalization  is 

5(0))  = a(o))E(o))  (2.3) 


and  the  frequency  dependent  conductivity  a is  a scalar  in  the 
one  dimensional  case. 

We  then  obtain  the  result 


AJ(t)  = Re{EQa  (o))e^‘^^}  (2.4) 

This  equation  relates  the  response  AJ(t)  to  the  excitation 
E^e^^^  directly,  and  the  conductivity  may  be  interpreted  as  a 
general  admittance  function  which  has  been  well  described  by  many 
researchers. In  particular,  Kubo^^  has  related  the 
general  admittance  function  to  the  response  function,  as  Fourier 
transforms  of  each  other.  Using  his  results  we  may  write 


AJ(t)  = a (t  - t")E(t')df' 
— 00 


(2.5) 
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0 is  the  response  function  and  eq.  (2.5)  says  that  the  sum  of 
excitations  that  act  on  the  system  from  the  present  time  back 
through  infinite  past  times  produces  a net  response  AJ(t),  the 
current  that  flows  in  the  plasma.  Kubo  derived  an  expression 
for  the  response  function  in  terms  of  an  equilibrium  ensemble 
average  of  a product  of  current  components,  namely 

3(t)  = ^ <AJ(o) AJ(t) (2.6) 


This  expression  assumes  the  plasma  to  be  a collisionless, 
single  component  system  of  volume  charge  density  n com- 
prised solely  of  electrons  having  charge  e.  The  angular 
bracket  represents  an  equilibrium  ensemble  average.  This 
is  an  extremely  important  point  in  the  interpretation  of  this 
formalism.  The  presence  of  the  external  electric  field  does  not 
push  the  initial  system  far  from  equilibrium  which  allows  us  to 
express  any  change  in  the  distribution  function  as  a linear 
functional  of  the  equilibrium  distribution  function.  This 
obviates  the  need  to  solve  some  difficult  transport  equation  for 
the  non-equilibrium  distribution  function,  and  allows  us  to  substitute  an 
easier  problem,  that  of  the  calculation  of  the  current  correla- 
tion function,  which  represents  the  effect  of  a perturbation  on 
the  originally  quiescent  system.  For  a three  dimensional  system 
the  frequency  dependent  conductivity  tensor  is  given  by  the 
Fourier  transform  of  eq.  (2.6), 


ne 

m*kT 


/ dte'^^'^^^^^  ^<P,,  (o)p  (t)  > 


o (uj) 

uv 


limit 

£-*■0 


(2.7) 
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We  have  recast  the  correlation  function  in  terms  of  the  charged 
particles  momentum  components  rather  than  the  current.  The 
exponential  damping  factor  is  a convergence  factor  that  eliminates 
any  persistent  oscillations  that  may  be  present  in  the  response 
function.  The  static  conductivity  is  the  zero  frequency  limit 
of  eq.  (2.7),  i.e. 


o (o)  limit  I dte“  ^^<p  (o)p  (t)  > 


(2.8) 


A further  comment  on  eqs.  (2.6)  through  (2.8)  is  required  to 
emphasize  that  these  expressions  result  from  a linear  approxi- 
mation to  the  Poisson  braclcet  of  the  distribution  function  and 
the  perturbation.  Retaining  higher  order,  non-linear  terms  would 
include  higher  inverse  powers  of  JcT  and  stronger  time  dependence. 
This  will  consequently  introduce  higher  inverse  powers  of  fre- 
quency in  the  frequency  dependent  conductivity  tensor.  Thus  the 
results  we  shall  present  in  this  paper  are  regarded  as  a linear 
approximation  which  may  require  further  refinement  in  the  future. 
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III.  THE  HAMILTONIAN 

We  Shall  consider  the  electrons  in  the  plasma  as  the  primary 
carriers  of  current,  and  shall  ignore  electron-electron  collisions, 
and  electron-ion  collisions.  The  Hamiltonian  for  such  a system 
of  N electrons  in  the  presence  of  an  external  line  charge  and 
current  source  may  be  written 


A(rj)  and  (ttr^)  are  the  vector  and  scalar  potentials  associated 
with  the  external  sources  and  the  electrons  in  the  plasma.  Let  us 
assume  there  is  a uniform  distribution  of  charge  in  the  z direction. 
Then  we  have  the  scalar  potential  expressed  as  a superposition  of 
plasma  sources  and  external  sources. 


N ^ 

<!>{]?)  = 'l>_(r)  + [ $.(r,  r .) 

° j=l  J J 


(3.2) 


where 


y <t>  (r)  = 
o 


_X_  6 (r) 
^o  ^ 


(3.3) 


and 


(r,0;rj,9^) 


, 6(r-r  ) 

— 1-  6 (9-0  . ) 

^o  ^ ^ 


(3.4) 


The  potentials  satisfy  Poisson's  equation  with  X and  e'  representing 
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the  charge  per  unit  length,  i.e.  dq/dz  for  the  external  source  and 
the  plasma,  respectively.  The  variables  r and  r^  represent  the  field 
and  source  radial  coordinates  respectively.  The  solutions  to  eqs. 


(3.2)  and  (3.3)  are 


$o(r)  = 


o 


(3.5) 


■ m(e-e.)]j  ,r>r. 


)j (r,0;rj,0j)  = 


(3.6) 


- ^ [irKr  /R)-  ^ ^(^j^'cos  [m(0-0  )]1 

o [_  m=l  ] ■’  J 


, r<r , 


where  R is  the  radius  of  some  conducting  shell  concentrically 
surrounding  the  line  charge  which  permits  the  potential  to  go  to 
zero  as  the  radial  variable  approaches  R.  R is  large  compared  to 
the  source  radial  dimension.  We  further  assume  that  there  are  n. 


electrons  having  source  coordinate  greater  than  the  field  coordinate 


(rj>r)  and  n2  with  source  coordinate  less  than  the  field  coordinate 


(rj<r).  The  total  scalar  potential  may  be  expressed  by 

ni 


1 r ^ 

4>  (r,  0;  r . , 0 . ) = - A - ■■  Xln(r/R)  + e'  ^ ln(r./R)  - 
3 J ^ j = l J 


(3.7) 


■ 2 » , r . 

+ n2e'ln(r/R)-  e'  [ [ - (^)  cos  [m(0-0.)] 

j =1  nT=l  3 
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The  vector  potential  is  treated  in  a manner  analogous  to 
the  scalar  potential.  It  satisfies  a vector  Poisson's  equation 
whose  solution  may  be  expressed  in  terms  of  Green's  dyadic. 
Presenting  just  the  results  we  have 


A(r,0)  = - 


^ n ev^ln 
ZTT  |_  O O 


(r/R)e  + 


<1 

1 °°  \ 
+ e'^  V jln(r  /R)  - I ~ if  )"’cos  [m(e-0  )J> 

j=l  J m=l  J ’ 


(3.8) 


+ e'ln(r/R)I^v.  {l-  ft  (J- cos  [ni(0-0.)]|] 

where  n^  represents  the  number  of  electrons  per  unit  length 
in  the  beam,  i.e.  dn/dz  and  v^  is  the  velocity  of  the  beam 
electrons.  The  coordinates  r^  and  0^  in  eqs.  (3.6)  and 
(3.7)  and  the  velocity  vector  v^  are  randomly  distributed  so  we 
shall  replace  them  by  their  expectation  values,  through  the  mean  fields 
approximations  (so  called  because  the  vector  and  scalar  potential 
fields  are  consequently  replaced  by  their  average  values) 


y f(v.)  = n.  / P (v  . ) f (v  . ) dv  ./  /P(v.)dv. 
j = l 3 j j J 1 J 


y fir.)  = n.  / P(r.)f(r.)dr./fP(r.)dr. 
j = 3 K 3 3 3 3 3 


(3.91 


"h 


y f(v  ,r.)  = nj^  / P(v  ,r  ) f (V  ,r  )dv  dr ./  / P(v  ,r  Jdv.dr 
j=]^  JJ  JJ  JJJJ  JJJJ 
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In  the  notation  of  eq.  (3.9)  P is  the  probability  distribution 
function  and  f may  be  any  arbitrary  function.  In  this  mean  field 
approximation,  the  average  vector  potential  momentum  <eA>  is 
assumed  small  compared  to  the  average  momentum  <p  > for  any 
electron  in  the  plasma.  This  zero  order  approximation  will  permit 
us  to  calculate  dynamical  quantity  averages  including  vector 
potential  effects  at  a later  stage.  The  initial  distribution 
function  we  shall  use  is 

P^(rj,Vj)  = exp  I -e  +e(t)^(rj)]}  (3.10) 

is  the  kinetic  energy  of  the  jth  electron,  e4)^  is  its  po- 
tential energy  with  respect  to  a shielded  line  charge  potential. 

For  a line  charge  of  radius  on  the  order  of  1 cm,  surrounded  by  a 
plasma  of  electrons  having  a Debye  wave  number  kj^  on  the  order  of 
10^  cm  ^ we  may  approximate  the  Debye  Huckl  potential  for  a system 
possessing  cylindrical  and  azimuthal  symmetry.  The  exact  solution 
is  derived  in  Appendix  A.  Quoting  the  result, 

o 

is  a modified  Bessel  function  and  is  related  to  the  ordinary 
Bessel  function  through  the  relation 

K^(x)  = J i (ix)  (3.12) 

Since  the  product  of  Debye  wave  number  and  radial  coordinate  for 
our  system  is  always  much  greater  than  unity  we  approximate  eq. 
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(3.12)  with  the  expression 


[r)  t ^ / e~^D^//r\ 

eo/8¥l^  V / 


(3.13) 


Calculating  the  average  vector  and  scalar  fields  using  eqs. 
(3.7),  (3.8),  (3.9),  (3.10)  and  (3.13)  requires  our  evaluating 

<ln(rj/R)>,  <(r/r^)"'>,  <(rj/r)"'>,  and  Because  the  Debye- 

Huckel  potential  of  eqs.  (3.11)  or  (3.13)  is  central  in  nature,  i.e. 


)o(r)  = <J>^(r) 


we  find 


<(r/rj)"'>  = <(rj/r)"*>  = 0 


(3.14) 


In  the  absence  of  the  vector  potential,  there  is  no  preferred 
direction  in  velocity  space  and  we  have 


<v>  = 0 

3 


I = (kjjR) 


(3.15) 


The  only  non-vanishing  average  is  the  logarithm  term. 

<ln(rj/R)>  = - (R-R^)"^  ^ (3.16) 


(3.17) 


6 = eX/e^/IrTkZ 

O D 


(3.18) 
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and  is  the  vanishingly  small  but  finite  radius  of  the  line 
source.  Substituting  the  results  of  eqs.  (3.14)  through  (3.16) 
into  eqs.  (3.7)  and  (3.8)  gives  the  average  vector  and  scalar 
potentials  that  shall  be  used  in  the  Hamiltonian  of  eq.  (3.1) 


- r 2 


U n V e j. . 

2 , o o o 1 , j V 

• + ^ ln(-:^)p  . + 

] TT  R ^zj 


ynve  2 -.r  ' 


- lAzl,  - n,(r  )Ne'(R-R„) 

0 j=l  '' 


(3.19) 


) + (R-Rq)  + 361 


The  number  of  electrons  n^^  and  that  enter  into  (3.19)  de- 
pends on  the  radial  coordinate  r ^ . Using  the  distribution  function 

of  eq,  (3.10)  and  the  fact  that  eiji  /)cT<<l,  we  find  that  the  electrons 

o 

for  all  practical  purposes  are  uniformly  distributed  throughout  the 
medium,  and  we  obtain 


n^(r.|  •v  np.LJR^  - r.^  ) 

n2(rj)  ■v  - R^^, 


(3.20) 


where  np  is  the  number  density  of  electrons  in  the  plasma,  and  is 
an  arbitrary  length  in  the  z direction.  The  average  number  of 
electrons  is  given  by 


<nj^(rj)>  = '^ri2(rj)>  np/2 


(3.21) 
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IV.  THE  LIOUVILLE  OPERATOR  AND  THE  TIME  PROPAGATOR 


The  Liouville  operator  for  a system  of  N electrons  is 


given  by 


L = - iUf?  ^ - 4 ] 

j=l  3pj  3rj  9rj  8pj 


(4.1) 


In  a cylindrical  coordinate  representation  with  azimuthal 
symmetry  this  becomes 


L - -i  I ^ 


3 r . 3p . 

1 jr 


(4.2) 


From  the  Hamiltonian  of  equation  (3.19)  we  immediately  find 


3H 

apjr 


(4.3) 


14  - + K2Pj2  + K3j  (i) 


(4.4) 


U 2n  2v  2^4 
O O o 


(4.5) 


(4.6) 


V - 

^3  * 2Tie. 


(4.7) 
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Substitution  of  eqs.  (4.3)  and  (4.4)  into  eq.  (4.2)  leads  to  an 
explicit  expression  for  the  Liouville  operator 


^ V:.  * ■'3]  r,  al-) 


(4.8) 


For  any  function  x of  momenta  and  coordinates,  p and  q,  which 
does  not  explicitly  depend  on  time, 


x[p(t),  q(t))  = exp(tiL)  X(p(tQ)  ^(t^) 


(4.9) 


"^12 

and  L depends  on  time  implicitly  through  p and  q.  It  should  be 
noted  that  L = L(p(t^),  q(t^)]  , that  is,  the  Liouville  operator 
is  expressed  in  terms  of  momenta  and  position  vector  at  some 
arbitrary  initial  moment,  t^.  The  operator  exp  (tiL)  is  Icnown  as 
a time  propagator  for  it  "propagates"  a function  at  some  initial 
arbitrary  instant  into  a later  time,  according  to  eq.  (4.9). 

Using  the  power  series  expansion 


exp  (tiL)  = 1 + tiL  + -^-(iL)  (iL)  + ... 


(4.10) 


allows  us  to  obtain  results,  repeating  the  application  of  iL  to  an 
arbitrary  function  of  a dynamical  variable  as  many  times  as  is 
necessary . 
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V.  THE  ELECTRIC  CONDUCTIVITY  TENSOR 


If  we  explicitly  utilize  the  propagator  of  eq.  (4.10)  in  the 


conductivity  tensor  of  eq.  (2.7)  we  find 


%v'“'  = ^ 

m kT  o 


X <exp(tiL) p^ (o) p^ (o) 


(5.1) 


Performing  the  indicated  operations  of  time  propagation,  ensemble 
averaging,  and  Fourier  transformation  we  find  the  components  to  be 


o („|  = ns!  r 2 . c.R 

rr  m i(jj  2 2m 


(1  - e-°2'») 

|(R  - 


2 - p 

= 2^  (4-)  r 

m 10)  I 


47ra2R^mkT 


(5.3) 


X I R^e"“2^  + 2Ra"^  (2e"“2*^  +1)  + 6a"^(e"°‘2^  _1)|  / 

/(R  - a“^  + a"^e”“2^)j 
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! 
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I 

i 

j 

1 
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o (UJ)  = !2® 

rz 


m kT 


[(2TTmkT)**  ^ |Ra‘^  (e"“2R  + i) 


(5.4) 


2a'^  (e*“2^_l)l  _ (i 

2 ; lu)  TT  \ 


a_R  -OL  R 

+ -4-  )e  ^ + 


1 “°‘2^ 
+ j (1  - e ^ 


/ ( ■*  - “2^  * “2^^  ) 


where 


a,  = y n V e 
1 o o o 


a_  = e(X  + n_e')/2Tre  RkT 
* 2 O 

and  R^,  the  radius  of  the  line  source  is  ignored  compared  to  R,  1 

the  radius  of  the  outer  conductor.  Examining  the  leading  terms  j 

of  the  conductivity  diagonal  elements  we  see  they  are  what  one 

would  expect  using  Drude  theory.  The  remaining  terms  involve  i 

higher  order  inverse  powers  of  frequency,  and  arise  from  the  * 

interaction  of  plasma  electrons  with  themselves  and  with  the  | 

external  line  source.  These  components  diverge  in  the  low 

frequency  limit  since  we  ignored  binary  collisions,  which  serve 

as  a damping  mechanism  for  the  system.  Thus  the  resistivity  of  ■ 

such  a system  in  an  applied  DC  electric  field  should  be  zero,  and  j 

its  conductivity,  infinite. 

The  z and  r components  of  current  that  flow  in  the  plasma 

are  given  by  a three  dimensional  extension  of  eq.  (2.5),  ■ 

' S 

3 

s 
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The  response  function  o (t)  is  found  from  the  Fourier  transform 

uv 

of  the  conductivity.  The  external  electric  field  may  be  written 
as  in  eq.  (2.2).  Performing  the  integration  in  eq.  (2.5)  yields 
a term  that  is  harmonic  in  time,  plus  higher  order  terms  that 
involve  a product  of  various  powers  of  time  with  exp(iwt).  Con- 
sequently the  current  increases  and  oscillates  in  time  and  will 

increase  indefinitely  without  approaching  a stationary  state. 

20 

This  result  differs  from  that  described  by  Bclt.scu  for  a homoge- 
nous plasma  in  an  externally  applied  field.  He  calculated  the  re- 
sponse current  from  a Green's  function  solution  to  the  Boltzmann 
equation  and  arrived  at  a stationary  state.  However  Balescu  has 
accounted  for  the  potential  energy  of  interaction  specifically 
through  an  iterative  perturbation  solution  for  the  Fourier 
transformed  Green's  function.  We  have  used  a mean  field 
approximation  to  arrive  at  a time  propagator  which  then  leads  to 
dynamical  quantities  that  monotonically  increase  with  time,  never 
reaching  a steady-state  value. 
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VI.  CONCLUSIONS 

The  correlation  function  approach  to  the  electric 
conduction  problem  of  a beam-plasma  system  is  successful  in 
predicting  the  frequency  dependent  conductivity  tensor  in  a 
mean- field  approximation.  The  local  behavior  of  the  conductivity 
is  ignored  in  this  calculation.  Such  an  effect  will  be  considered 
in  a future  paper.  The  effects  of  electron-electron,  electron-ion,  and 
electron-neutral  binary  collisions  have  also  been  ignored  and  it  is 
hoped  that  this  too  shall  be  studied  at  a later  date.  The 
general  behavior  of  the  conductivity  components  follows  the 
classical  Drude  theory^ results  with  the  interaction  between 
beam  and  plasma  showing  up  as  higher  order  terms  in  inverse 
powers  of  frequency.  The  off-diagonal  component  does  not 
possess  Drude-like  behavior  as  its  leading  term,  however. 

Thus  the  response  current  in  the  plasma  increases  with  time 
indefinitely,  and  a steady-state  or  purely  oscillating  current 
is  not  obtained  within  this  model. 


] 

i 

i 

I 

I 

I 
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APPENDIX  A.  THE  SHIELDED  LINE  SOURCE  POTENTIAL 


According  to  the  Debye-Hiickel  theory  as  presented  by 
Jackson^  we  imagine  the  line  charge  surrounded  by  a specified 
distribution  of  electrons  with  a uniform  background  of  positive 
charge.  The  distribution  function  for  the  electrons  will  be  the 
canonical  one  described  in  eq.  (3.10).  Then  Poisson's  equation 
is  written 

r\  ^ _ 

-6e 


= A 


r)  "o^  r 
^ ■ T-  L' 


"b-  1 


(A.l) 


The  first  terms  on  the  right  hand  side  is  the  charge  distribu- 
tion of  the  external  source,  the  second  represents  the  electron 
plus  uniform  ion  distributions.  6 is  just  the  Boltzmann  factor 
(kT)  For  a plasma  having  a temperature  of  a few  thousand  degrees 

with  a Debye  wave  number  of 

kp  lO^cm  ^ , the  condition 
etfijr)  /kT  < < 1 

is  easily  met.  So  we  may  linearize  eq.  (A.l) 


, d<()  2 , 

d?  (^dF^  - ^’^D  <^0=  - — 

o 


(A. 2] 


where 


(A. 3) 

is  the  Debye  wave  number.  We  have  assumed  cylindrical  symmetry  to 


arrive  at  eq.  (A, 2).  Its  solution  is  found  in  Morse  and  Feshbach. 

O 
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(A. 4) 


NSWC/WOL  TR  77-146 


APPENDIX  B.  AVERAGE  VALUE  OF  THE 
LOGARITHMIC  TERM  IN  THE  MEAN  FIELD  APPROXIMATION 


The  average  value  of  the  logarithmic  term  in  the  expressions 
for  the  scalar  and  vector  potential  of  eqs.  (3.7)  and  (3.8)  is 
found  from  eqs.  (3.9)  and  (3.13). 


:ln(rj/R)>  = / exp  | -g  j^6exp  (-jc^r ^ ) /r  x 

^o 

c ln(rj/R)drj  / / exp  | -6  pexp  (-Ic^^r  ^ ) /r  J dr^ 


(B.l) 


The  limits  of  integration  reflect  the  upper  and  lower  bounds  on 
the  radial  coordinate,  R representing  the  radius  of  the  outer  con- 
ducting shell,  and  R^  the  small,  but  finite  beam  radius. 

Inspection  of  the  exponential  argument  reveals  that  a lineari- 
zation is  a suitable  approximation  for  the  typical  beam-plasma 
systems  that  are  of  interest  to  us,  i.e. 


MAX|e())Q(rj)  1 /IcT  <<  1 

Retaining  the  linear  first  order  term  in  the  expansion  of  the 
numerator  and  unity  for  the  denominator  in  eq.  (B.l)  yields 

<ln(rj/R)>  =»  - (R-R^)"^|^^ln(I^/R)  +(R-Rq)  + e63;|  (B.2) 

where 

I = /exp{-  )c  r.}  ln(r./R)/r.’  dr.  (B.3) 

Rq  3 3 3 3 

The  integral  in  eq.  (B.3)  may  be  expressed  as  a combination  of  in- 
complete gamma  functions  if  we  use  the  leading  term  in  a series 


NSWC/WOL  TR  77-146 
approximation  for  the  logarithm.  Then 

lnx*x-l,  0<xl  2 
We  find 

I = - (k^R)'S(^,kjjR) 

- R**  - (kj^R)’^*  Y(j,kj5R^)J  (B.4) 

The  incomplete  gamma  function,  y is  defined  in  terms  of  the  ordinary 
gamma  function  r(c)  and  the  complementary  incomplete  gamma  function 

r(^,5) 

r(0  = r(c,c)  + Y(c,C)  (B.5) 

r(c,C)  has  an  asymptotic  representation  (since  k^^r  >>  1) 


i 

(B.6)  1 


i 

i 

I 


Conseauently  we  obtain 
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